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I. INTRODUCTION 

The present paper continues work started in Refs. where a boson-expansion tech- 

nique, which is widely used in nuclear structure theory, was applied to models of quantum 
field theory (in particular, the linear sigma model |3|). This technique appears to be a 
promising tool for building symmetry-conserving and non-perturbative approaches. It was 
demonstrated that the Holstein-Primakoff mapping (HPM) is able to systematically classify 
the dynamics according to the 1/N expansion. For example, a tttt scattering equation to 
lowest order in 1/N (which is equivalent to the random phase approximation (RPA)) was 
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evaluated in Ref. [3] and fulfills all required properties - unitarity, the Ward identity and 
standard RPA sum rules. 

However, the perturbative boson expansion in general and the HPM in particular rely 
on the bozonization of pairs of particles. Thus, in |2j single-pion images were missing 
after the mapping. Their lack proved to be an obstacle for defining unambiguously the 
two-point Green's function for the Goldstone mode. Another problem is the incorporation 
of single-particle states in the HPM which is called for in finite-temperature applications 
of the approach. Here, the absence of single-particle modes causes serious difficulties in 
the treatment of thermal ensembles since the single-boson (fermion) density of states plays 
an important role. These shortcomings were overcome in Ref. j^J where the bosonic HPM, 
originally introduced in p, was extended to accommodate a single-boson mapping. When 
applied to the bosonic O(N) anharmonic oscillator to leading order in the 1/iV-expansion, 
the extended HPM (EHPM) was shown to render N uncoupled Goldstone- as well as RPA- 

n 

phonon modes. Subsequently it was found in Ref. p that, without single-boson mapping, 
the thermal Goldstone modes were missing from the spectrum of the O(N) model to leading 
order in 1/N. 

In spite of the obvious success of the formalism its full power should be revealed by 
working out the next-to-leading order since this allows to evaluate n-point functions unam- 
biguously. This is the main aim of the present work. Although the 0(N + 1) anharmonic 
oscillator is a purely quantum mechanical system its treatment in the 1/N expansion, which 
is largely analytical, provides valuable insight for the use of non-perturbative methods in 
quantum field theory (QFT). 

The paper is organized as follows. Based on the Hamiltonian of the 0(N + 1) anhar- 
monic oscillator, its perturbation series in powers of 1/yN up to next-to-leading order is 
constructed in Sect. 2 within the EHPM framework. From its diagonalized form the validity 
of the Goldstone theorem is inferred and the correction to the condensate shift is computed. 
The validity of the Ward identity in the same order is proved in Sect. 3. A summary and 
conclusions are given in Sect. 4. 
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II. PERTURB ATIVE EXPANSION OF THE HAMILTONIAN TO ORDER 1/N 

The properly scaled Hamiltonian of the anharmonic oscillator with 0(N + 1) symmetry 
broken down the O(N) group reads: 



p2 p2 2 

H= -^ + ^ + — 
2 2 2 
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- VN V X a , (1) 



where we have chosen a notation reminiscent of the linear a model ^3], to emphasize the 
analogy to QFT. The variables X n and X a are defined as: 

1 1 N 1 

X n = (a + a + ) = + 4) , X a =^=(b + b+), (2) 

and P n and P a denote their respective conjugate momenta. The frequency, £ a , of the mode 
X a will be defined later. 

Diagonalization of the Hamiltonian leads to the excitation energies of the system. A 
standard approximation scheme in many-body theory consists of splitting H into a mean- 
field part and a residual interaction. The mean field is obtained selfconsistently in the 
Hartree-Fock (HF) or Hartree-Fock-Bogoliubov (HFB)approximation which becomes exact 
in the large- iV limit. It was shown in Refs. y], 0, y] that, for any finite N, the HFB 
approximation leads to a non- vanishing 'pion mass' in the exact symmetry limit (rj = 0). 
Hence, there are no massless Goldstone bosons in the mass spectrum. However, in the 
large- N limit, the Goldstone mode is recovered. This stems from the fact that the Fock- 
term is subleading and therefore vanishes. The Goldstone mode is provided solely by the 



Hartree approximation. Thus, the state o 



treatment of the residual interaction. In 



the 'Hartree pion is suitable for a perturbative 



^], such a state was introduced in the framework 
of the EHPM and allows to preserve the original single-pion degrees of freedom. The EHPM 
method consists of the following mapping: 



(aa)i = \/2N + 4(n + m)A , 
(a + a)j = 2n + m , {a + a + )i = (aa)/ , 
(a,)/ = -\/2N + A(n + m)T N (m)ai + 2aiAT N (m) , 
(af)! = ■ (3) 



1 We will use Goldstone mode and 'pion' interchangeably 
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Here N is an integer (the number of Goldstone modes), n = A + A, m = Yli a t a i an d is 
given by 

T N (m) 



m + N-2 ] 1/2 

(4) 



_2(2m + N)(2m + N -2) A 
Thus, instead of the original pion aj, we have an ideal boson at which is a Hartree pion |4|. 

A vacuum state of the system is a coherent state |4] which incorporates spontaneous 
symmetry breaking via the condensates for the A and b bosons 

=exp[(A)A + + (b)b + ]\0) . (5) 

The mode a, on the other hand, is not allowed to condense. The state is annihilated by 
the shifted operators A and f3 

A = A-(A), (3 = b-(b), (6) 
where the condensates (A) and (b) are scaled with iV according to 



(A) = d ) j^, (b)=s^. (7) 

After a redefinition of A and (3 a systematic 1/N expansion of the Hamiltonian can be 
performed 

h = nh q + Vnh 1 + h 2 + ^h 3 + —h 4 + --- . (8) 

VN N 

where H denotes the mean-field Hamiltonian. The fact that iV is arbitrary insures the 
preservation of the symmetry of H at each order. The dynamics up to order iV -1 is specified 
by the Hi terms in (jHJ) with i — 1, ... ,4 which are given explicitly in Appendix A. 

The coherent ground state is determined by a requirement that the ground-state 
energy NE = (^f\H\^f)/ (^|^) is minimized (note that only H contributes this expectation 
value as discussed in Appendix A). The minimization procedure with respect to s and d then 
gives the following two coupled gap equations: 

UJ s 



2udVT+d^ + A(Vl+rf2 + d f = o } (9) 



where the gap parameter A is given by 



UJ UJ 2 
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In fact, this is equivalent to the requirement that Hy — 0. 

With the gap equation 0, and the masses of the Hartree pion (£ n ) as well as the sigma 
boson (S a ) (see Appendix A for details) 

S T = u +4gs + — = - , 

C-7T S 

Si = u 2 + 8gs 2 + ^- , (10) 

the Hartree ground-state energy and the quadratic part H 2 of the Hamiltonian expansion (jHJ 
can be written in a more compact form which is particularly suitable for the further discus- 
sions 



S 2 

TT 



5 2 

s~~ 3s 



Si 



2 1 

s + s~ 



4£ 



(id 



2 ' 



H 2 = S W J2 a t^ + S a (3 + (3 + 2S n B + B + ^(B + + Bf 



+-^=(P + + P)(B+ + B) + ^. (12) 



The operators B , i3 appearing in ()12j) are related to A and ^4 a through unitary transfor- 
mation (see (|A5|) - (|A6|) in Appendix A). 

The gap equation and the mean-field masses of the Hartree pion and the sigma- 
boson (fTUj) can be presented diagrammatically as shown in Fig. ^1. 

From Eq. (fT2*|) it follows that H2 describes several modes of excitation above the Hartree 
ground-state energy. On the one hand, one has independent Hartree pions which have the 
Goldstone character, i.e. their frequency S n vanishes in the exact symmetry limit (rj = 
0), and for a finite condensate (s 7^ 0). Besides the Goldstone modes there also exist 
other excitations which can be made explicit by diagonalizing the remaining part of H2 in 
RPA 

The next two terms, H 3 and H4, in the Hamiltonian expansion (jHJ) renormalize the energy 
of the modes obtained after H 2 diagonalization and represent anharmonic corrections. The 
c-numbers in H 2 and H4 contribute to the ground-state energy. 

To diagonalize H 2 , Hz and we use the method described in detail in Ref. p|. We 
have used it previously jSj to diagonalize next-to-leading order terms in the 1/N expansion 
of the Hamiltonian of the Lipkin model jj. In accordance with the prescriptions of Ref. []| 
the diagonalization of (jHJ) within perturbation theory can be carried out in two steps. At 
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FIG. 1: a) Diagrammatic representation of the mean- field results for pion and sigma bosons. The 
solid line denotes the self-consistent pion propagator, the wavy line - the sigma propagator. The 
dashed (dotted) line is the two-point Green's function of the bare ir (a), b) The sigma-boson 
propagator in the RPA order (a thick wavy line). 

the first stage RPA bosons are introduced by means of a linear canonical transformation. 
This diagonalizes the sigma part of H^. These RPA bosons are then used as a basis for a 
systematic diagonalization of the higher-order terms H 3 and H4 and a remaining part of the 
Hamiltonian should be written in terms of RPA bosons. The further diagonalization can be 
accomplished by means of a unitary transformation which removes all off-diagonal terms in 
any order leaving the Hamiltonian finally as a function of a new " diagonal" bosons. 

Following this procedure we define first a new RPA boson operator and Q v via a 
general Bogoliubov rotation which mixes the sigma boson, /3, and pion pairs £>: 

Q+ = x v p + + y v p + UuB + +v v B. (13) 

By the equation-of-motion method we derive RPA equations for the amplitudes 
x u , y v , Uu, v u . The corresponding system of equations can be reduced to the following 
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two equations for the RPA frequencies Q v [y = 1,2) (see also Appendix B): 

4<? 



n 2 + n 2 = 4% + % + 



(14) 



Formally, the solution of (|14|) can be written as 



2 + 

S 



8gs* 



i _ ?y 



(15) 



The inverse transformation is given in Appendix B. 

Next we need in a sigma-boson propagator in RPA, to order, i.e., 1/N. It can be evaluated 
by solving the Dyson equation with the mass operator taking into account pion fluctuations. 
The diagrammatic representation is shown in Fig. [I]p. The corresponding Dyson equation 
has the form 



D a {k 7 



+ 



k 2 - E 2 k 2 -E 2 



32g 2 s 2 )E 7T7T (k 2 )D a (k 2 ) 



(16) 



where S nv (k 2 ) is an infinite sum of two-pion loop diagrams 



^TTTT {k 



dp 



27T (p* - E$){\p - k\ 2 - e$ 



-Aig) 



dp 



2tt ( P 2 - S 2 )([p - k} 2 - E 2 ) 



+ 



E^(k 2 -AEl 



i-4ig) 



[E„(k 2 -4E; 



+ 



(17) 



E^k 2 -4El-4g/E„) 

The additional factors 1/N at each vertex are precisely canceled by the factor iV appearing 
due to summing over N internal pions in each bubble. Inserting the explicit form of E wlT {k 2 ) 
into Eq. (fTHjl one obtains 



D a {k 2 



(18) 



k 2 -E 2 + 2,2ig 2 s 2 E^(k 2 ) 

Now the RPA frequencies Q u can be obtained by solving [D a (Q' 2 ,)}^ 1 = 0. They are identical 
to those in Eq. ()15j) which were obtained by using the equation-of-motion method. In terms 
of 0^2 the sigma-boson propagator reads 

. k 2 - 4E 2 - 4g/E n 



D a (k A 



(k 2 -ni)(k 2 -ni 



(19) 
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It should be mentioned that in this approximation the condensate s does not receive any 
contributions from RPA fluctuations. Moreover, if k 2 —>■ Q 2 , the sigma-boson propagator 
becomes 

D - {e) ~ ¥%i • (20) 

where the field renormalization constants Z v are given by 

3, = - j&)/(n£ - n£) = . (21) 

Here and hereafter the notation v' means that v' ^ v. 

Now, the H 2 part of the Hamiltonian can be written in diagonal form in terms of Hartree 
pions and RPA bosons 

H 2 = S^ a t®i + ^QtQv + Ql + Q2 - £ w . (22) 

i=l i=l,2 

The last two (c-number) terms in (|2*2*|) represent the RPA corrections to the Hartree ground- 
state energy NE . 

Having diagonalized H 2 we now turn to if 3 and H4 and express them in terms of the 
RPA phonons. The explicit expressions are given in Appendix B. To diagonalize we use the 
q _ analog of the classical Bi rkM -G_ approa^ Q Q as was done in [8] . The 
diagonal terms in i? 4 (Appendix B) provide anharmonic corrections of order iV -1 to RPA. 
On the other hand, all terms in H% include only odd powers of Q u , and have, therefore, 
vanishing diagonal matrix elements in the RPA basis. They first contribute in second-order 
perturbation theory and are thus of the same order as the diagonal contributions of H4. 
Therefore, these parts of the Hamiltonian must be considered together. 

As has been mentioned, it is possible to define a perturbative, unitary transformation 
which removes all off-diagonal terms in any given order of the l/N- expansion. At the order 
of interest, this transformation has to be chosen such as to eliminate the if 3 part of the 
Hamiltonian expansion (jHJ), thereby shifting their effects to renormalized diagonal terms of 
F 4 . 

To this end, we define new "diagonal" pion and sigma operators (ctf)d an d (Qt)d which 
are obtained via a unitary transformation from the «j and Q u operators such that 



[ol. 



-) d = e^^afe^ , (Q+) d = e-^Q+e^ . (23) 
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The Hamiltonian is then rewritten in terms of these "diagonal" operators 2 as: 

H = e -is d /VN Hde i Sd /VN = Hd + -L[ H ,iS] d + ±-[[H,iS],iS\ d + ■■■ . 

v N ^-iv 

Using the Hamiltonian expansion (JHJ), we obtain after some algebra, in the order of interest 

H = NH + (H 2 ) d + J= {[H 2 , iS] + H 3 } d 

v N 



+J-{[[H 2 ,iS],iS} + 2[H 3 ,iS} + 2H i } d + ... . (24) 



Therefore, the condition for eliminating the off-diagonal terms of order of 1/vN can be 
easily read off as 

[H 2 ,iS] d = -(H 3 ) d . (25) 



and completely determines the operator iS d . Using this result in Eq. (j2Tj) we obtain after 
some algebra (Appendix C) the full diagonal part of the sigma-model Hamiltonian in next- 
to- leading order of the 1/N expansion: 



H = N(H + Ql + ^ N ^ ) + (^ + 1^) J>^) d 

i 

+ + -^ Q '»)(QtQ»)d + Jj{\j2 Vnr-^at afoot 

v i,j 

+ J2 V™»->™MQtQv<*i + \Y v ^ v ^ v/ QtQ + u>Qv>Qv} d ■ (26) 

i,v u,v' 

Besides RPA corrections to the Hartree ground-state energy it contains all anharmonic 
corrections caused by residual interactions between different modes as well as higher-order 
contributions £' n , Q! v to the masses of pion- and sigma-bosons, respectively As can be seen, 
there are three sources of residual interaction: the pion-pion term ~ V^^,^, the pion-sigma 
term ~ K-ov->7rovj an d the sigma-sigma term ~ V auCT ,^ ava ,■ The explicit form of these 
couplings as well as of S' n and Q' v are given in Appendix C. It is evident from (jC7|) and (jC6j) 
that the 1/N order correction £' n to the pion mass and pion-pion interaction term V^^^ 
vanish in the exact symmetry limit rj = 0. This implies that the pion keeps its Goldstone 
character in next-to-leading order of the 1/N expansion. 
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Here we have adopted the following notation. For any operator F, depending, e.g., on O and + , Fd is 
obtained by the replacement O — > Od and + — > Oj~, i.e. F — e lSd Fde~ lSd . Of course, one has e lS — e lSd . 
We also use the following appreciation for a product of operators: adbd ■ ■ ■ Zd = {ab . . . z) d - 
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Next-to-leading order corrections also contribute to the sigma-boson shift (|7j). The cor- 
responding correction reads 

where (...) denotes the expectation value with respect to the vacuum of "diagonal" bosons. 
Expressing + (3) through these bosons we obtain 

2< 8gs / 3 3 \ 

III. WARD IDENTITY IN NEXT-TO-LEADING ORDER 

Aside from the Goldstone character of the pion the well-known Ward identity relating the 
pion propagator to the sigma-boson shift [13J should also remain valid. In this Section we 
prove that the Ward identity is indeed fulfilled in next-to-leading order of the 1/N expansion. 

The Ward identity means that the pion propagator with 1/N corrections to the pion 
self-energy and the sigma-boson shift value calculated in the same order should satisfy the 
following relation: 

M" 1 = jrh? ■ (28) 

To prove Eq. (}2*%|) we will use a diagrammatic technique to obtain the pion propagator in 
the order of interest. A necessary step is to calculate the 7T7r-scattering T-matrix which 
contributes to the pion self-energy at order 1/N. There exist two types of ttti- interactions: 
those involving contact interactions (first two diagrams in Fig. EJ) and those arizing from 
sigma-boson exchange (remaining diagrams). It is important to stress that, in the 1/N 
expansion, diagrams of the same topological structure can contribute in different orders. If 
% 7^ j in Fig. El then there is no summation over an internal pion in the s- and w-channels. 
Adding a new pion bubble increases the order of the diagram by 1/N. Therefore, in this 
case an interaction exists only in t-channel to leading order. On the other hand, in the case 
i = j all the three channels s, t, u contribute to the scattering process. 
The first two diagrams give the following contribution to the T-matrix: 

T cont (/ 2 ) = 8g {1 + {-Aig)E^{l 2 )} = 8g [2 - ^ 2 _ * , (29) 
where I = k\ — A given channel is specified by a fixed value of I 2 : t-channel I 2 = 0, s-, 
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• P2 Pi 




+ 



+ 



u and s 
channels 
if i = j 



FIG. 2: Diagrammatic representation of tt-tt scattering. The wavy line denotes the sigma boson 
RPA propagator while U nn contains an infinite sum of interacting pion loops. 

u-channel I 2 = AS 2 . From the a— exchange diagrams we have 



T exch (/ 2 ) = -56^ 2 ^ CT (Z 2 ){l + (-4^)^(Z 2 )}' 

56g 2 s 2 (l 2 -AS 2 ) 2 

(I 2 - AS 2 - Ag/S^)(l 2 - £l\){l 2 - Of.) 



Thus the full T matrix in a given channel reads 



T^i} 2 ) — T cont (/ 2 ) +T exch (7 2 ) — 8g 



(I 2 - AS 2 )(l 2 - S 2 ) 

(i 2 -n 2 )(i 2 -n 2 ) 



(30) 



(31) 



As was already mentioned, if i ^ j, only the interaction in the t-channel should be taken 
into account, hence I 2 = 0. Then 

32^ 



T 



7^(0) 



On the other hand, if i — j, 



T-KiiKi 2!Z^ r7r (0) -\- T nn (AS n j 



64gg 



Dividing the T matrix by a factor 4£^(l+<5y) which arises from the relativistic normalization 
condition, we obtain the matrix element of the pion-pion interaction V^^^ (jC6|) calculated 
in the previous section in the framework of the EHPM. 

We are now able to calculate the pion propagator to order 1/N. The diagrammatic 
representation of Dyson equation and the condensate is shown in Fig. El 
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FIG. 3: Upper part: Dyson equation for the pion propagator (doubled line) to order 1/N. Lower 
part: Dyson equation for the sigma-boson condensate value in the order 1/N. 

As one can see there appear two additional diagrams as compared to Fig. involving a 
sigma bubble and T^. Thus, to order 1/N, we obtain two coupled equations for the pion 
propagator D n (k 2 ) and the sigma-boson shift s + N~ 1 s' 

[iD^k 2 )]- 1 = -k 2 + u 2 + Ag{s + N~ l s') 2 + ^ J ^{p 2 ) 



(32) 



{s + N-\s') 



s + N^'s 



dp 



2tt 



D^p 2 



4<? 



Nu 2 



(s + N-'s 



dp 
27 



D a (p 2 ){3 + 2(-4ig)Z 7T7T (p 2 )} 



-%(s + N~hr 



(33) 



where D a (p 2 ) is the sigma-boson propagator in the RPA order. Except for terms involving 
the pion propagator D^p 2 ) all other pieces in (J3*2)l and (J3l?j) can be calculated analytically. 
After some algebra we obtain 



V 



s + JV-V 



lu 2 + Ag(s + N- 1 s , y + 4g 



dp 



2tt 



D«(p 2 



2 A 

n [n^nj-n 2 ) n 2 (n,l-n 2 )_ 



Ag_ 
' N 



< - n? 



(34) 



This proves the Ward identity (|28|) in next-to-leading order in 1/N expansion. 

We show that the condensate s' in Eq. (J3ll)) are equal to that in Eq. (|27|) which is calculated 
with the EHPM technique. To do so we should express the pion propagator (J3*2j) in terms 
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of the pion mass in next-to leading order, (E^ + iV 1 S' n ) , which is the solution of 

[DAk 2 )V = • 
In terms of this mass, the pion propagator reads 
[^(A; 2 )]- 1 = -k 2 + (S n + N~ 1 £' 7T ) 2 

-| / [TU\£,-P?)-TU[k-Pm ■ (35) 



Putting k 2 = we obtain 



[^(O)]- 1 = (S^ + N-'S!^ 
8g f dp i 

= (£ w + N-'Sl) 2 



N 



C-7T 



_Q 2 Q 2 tt^Ql-tt 2 ) tt 2 {tt 2 - tt 2 ) _ 



(36) 



Using this result in the lhs. of Eq. (|28j). one recovers the result in Eq. (J2ZJ). 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have applied the extended Holstein-Primakoff boson mapping to 
the bosonic O(N) anharmonic oscillator which mimics some essential features of the linear 
sigma model of QFT. The approach that we use fundamentally relies on the bozonization of 
the generators of the Sp 2 <8> Ni semidirect product group and is very appealing from many 
aspects and particularly from its symmetry-conserving character. The extended boson HP 
mapping exploits both single-boson and two-boson images. The former seems to be of crucial 
importance to correctly embed the correct quantum statistics in the mapping [Q| . Moreover, 
it permits unambiguously to define the two-point function for the Goldstone mode. 

In Refs. 0,0], the O(N) anharmonic oscillator has been treated within the EHPM for- 
malism to leading-order in the 1/N expansion which is equivalent to the Hartree + RPA 
scheme. In that approximation the system appears as a system of iV Goldstone modes (pi- 
ons) and two RPA modes which are superposition of the initial sigma-boson and auxiliary 
"two-pion" modes. 
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In next-to-leading order the above excitations are coupled with each other. We have 
calculated the corresponding coupling matrix elements 3 . Furthermore, next-to-leading order 
corrections to pion- and sigma- masses as well as to the condensate shift were obtained 
through perturbative diagonalization of the Hamiltonian using single-pion and RPA bosons 
as a basis. We have proved that in next-to-leading order the Goldstone theorem remains 
valid and the Ward identity is satisfied. These results constitute an encouraging step towards 
the application of the EHPM techniques to field theoretical modes, such as the linear sigma- 
model. 
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APPENDIX A 



In this Appendix the boson expansion series up to next-to- leading order in 1/N for the 
Hamiltonian (1) is presented. 

In terms of the shifted boson operators (3 and A the X 2 and X 2 take the following form: 

Xl = Ns 2 + VNs^-(P + + 0) + ^-(/? + + (3) 2 , 

X~l = ^-\N(l + 2d 2 )/2 + V2Nd(A + + A) + 2A + A + m 
2uj I 



+V2NK(A + + d^/N/2) x 

2m + 2A + A + V2Nd(A + + A) 
+ N^K 



1/2 



h.C. 



(Al) 



where K — 1 + d 2 . After a Taylor expansion of a square root in (jAlj) up to order l , the 



3 For the present, purely quantum-mechanical model, these matrix elements play a role of four-point Green's 
functions for QFT in four space-time dimensions. 
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operator X 2 is given by 



Xl = N^-(VK + d) 2 + N 1 / 2 iVK t^} (A + + A) 
4w v 2luV2K 



(A + + A) 2 - - (A + - A) 2 + 4(m - 1) 
K 



K + d 



8uj 



(A + + A) 2 - - (A + - A) 2 (VK - d) 2 
K 



+N- 1 ' 2 



(A + + A) 
d 



(A + + A) 2 \- - (A + - A) 2 + 4(m - i; 
K 



64^3/2 



CO 



5(A + + Af\- - (A + - A) 2 + 4(m - 1) 
K 



1 



(A + + A) 2 - - (A + - A) 2 + 4(m - i; 
K 



4(^ 2 -i; 

K 



+ h.c. 



(A2) 



Substituting this expression into Eq. (JSJ) leads to the required 1/N expansion of the model 
Hamiltonian. The first term NHq yields the ground-state energy in the Hartree approxima- 
tion 

(V\H\9) 



En 



CO 



N{\V\W) 



-(1 + 20 + 



+9 



u 2 s 2 



K + d) 1 
~4u7 2 



+ s 4 + 



K + dy 
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(A3) 



Minimizing E determines the values of d and s as the roots of the coupled gap equations Q. 
The next term H\ is linear with respect to A and f3: 

{A + + A) 



Hi 



V2K 



2(vdVK + A(VK + d) 2 
V 



(A4) 



and vanishes at the minimum of Eq. 

For the following it is particularly suitable to define a new operators B + and B 

A + = G+B + + G-B . 



(A5) 



In order for the B + , B operators to obey the same algebra as A + , A the coefficients G + and 
G- have to be chosen such that G 2 + — G 2 _ = 1. Explicitly, they are given by 

1 



r— 1 



(A6) 
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and we have 



(A + + A) = (B + + B)Vk , (A + -A) = (B + - B)/VK 

In terms of B + and B the quadratic part of H now reads as 

+ d 



(A7) 



Ho 



uj + A 



K 



+ 
+ 



uj 2 - El ujA 2gs 2 



4£cr 2S a S a 
uj 2 + £ 2 a , ujA 2gs 2 



(2(3 + (3 + l) 



+ ^K + d y {B+ + B)2+ 2. B+B 



2uj 2 

2gs(VK + d) 2 



K + d) 2 
(f3 + + (3)(B+ + B) . 



From this expression, and more precisely from the coefficient in front of the term £\ a^ai, 
one can deduce the existence of N uncoupled modes that are nothing but the Hartree pions 
with mass S n given by 

^K + d uj 



£ n = uj + A 



K + df 



El = uj{uj + 2A) = uj 2 + Ags 2 + j- . 



(A8) 



Furthermore, demanding that the bilinear part of H2 in the f3 operators is diagonal, i.e. 

uj 2 - Si ujA 2gs 2 



4£„ 



we obtain the sigma-boson mass £ a 



2£ n 







Si = uj 2 + 12gs 2 + ^-=Sl + 8gs 2 . 



(A9) 



With these definitions of S n and £ a the quadratic part H 2 of the Hamiltonian can be recast 
in the compact form given in Eq. f)12j) . 
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The next-to-leading order terms H 3 and H 4 in the Hamiltonian expansion look as follows: 

x (4(771 - 1) + (8+ + /3) 2 - (8+ - Bf\- 

K 



(A10) 



x {^(/3 + + /3) 2 - -^-(/3+ - Bf + + /3) 2 } + h.c. 



//; = ^(U + T-)('» "I) 2 



+-|^(B+ + B)(;3++;3) +! ,f-!- + i 



X 



^(/3+ + /3)(/3 + + /3) + r 7 (l + i') 



32K 
8gs 



x{(B + + B) 2 -i(B+-B) 2 } 



^ 2 - 1 Uh.c. (All) 



iK 2 £^u 8K 2 

APPENDIX B 

In this Appendix some details concerning the Bogoliubov rotation (|13j) are presented. 
We start with the following definition of the RPA excitation operator: 

Q+ = x u {3+ + y u (3 + u u B + + v v B . 

Using the equation of motion 

(RPA|[5Q„, [tf 2 ,Q+]]|RPA> = ^(RPA|[5Q„ Q+]|RPA> , (Bl) 
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where the RPA ground-state is determined as usual by the requirement that Q^|RPA) = 0, 
we obtain the following characteristic equation for RPA frequencies: 

# - nl [a£1 + El + + < (El + j- ) = . (B2) 

The solution is given by 

s 1 4<? 



Furthermore, using the relations in Eq. (fTl| it is easy to verify that the RPA frequencies 
obey the following condition: 

(nl-s^-nl)^ 3 ^, (B3) 

and we shall assume that flf > 

Using the orthonormalization condition of RPA states (RPA|Q„Qj|RPA) = 5 uu i , one 
can find the RPA amplitudes x u , y u , u u , v v and express (3 and B operators through Q v . 
However, as is seen from Eqs. (jAlOj) and (jAll)) it is more convenient to use the inverse 
transformation Q u —>■ j3, B in the following form: 



2E 2£ 
(B + +B) = J^ ( p 1 (Q+ + Q 1 ) + ^l^ 2 (Q+ + Q 2 ), 

{B + -B) = J^ ( p 1 (Qt-Q 1 ) + x [^MQt-Q2), 



(P + + P) 



■][^MQt + Qi) + \[^MQt + <h), 




{(3 + -(3) = - A / T ±0 2 (Q+-Q 1 ) + A /_f0 1 (Q+-Q 2 ), (B4) 



where 01 = v 02 = v ■ 

Now we are able to write down the E 3 and E^ parts of the Hamiltonian in terms of RPA 
bosons Q u and Q+. For the E 3 part one obtains 

E 3 = H™ + E™ , (B5) 

with 

HT = ^ {^(Of + ft) + + ft)} (m - 1) , (B6) 
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= J2 A AQt + Qvf + a v2 {q + v + g„) 2 (g+ + q,) 

v=\,1 

+A u3 (Qt - Qu) 2 (Qt + Qv) + A vA {Qt - Q V )\Q+ V , + Q v >) 

+A v ,{Qt - Q v ){Qt + Qu)(Qt> + Qv) + (h-c.) , (B7) 

where coefficients A vi (i — 1, ... ,5) are 

_ <f>v ( <^K p 



A 



v3 



/ R v K 




and 



/07 2 <^ / ^ 



4 V u) n 7 4 

Note, that H™ is proportional to (m — 1) (m - is the pion number operator) whereas ifg' 

only includes the RPA boson operators Q u , Q+. 

Finally, the if 4 part of the Hamiltonian expansion reads: 



with 



H A = H? + H? + H% a (B8) 



1 (a ,2 , 9 \ , -, 



HT = x[£«d < + ^-)(m-l) \ (B9) 



x{h^) + ^m} {2QiQ ' +1){m - 1) 

+ (off - diag.) , (BIO) 

H ™ = ^{tk (jt + ^k) Tt + tk (1: ~ ^k) k +Uu ] 

x(QtQtQuQu-2QtQu) 

+ 2 \AK V ^2 2£*K ) 1 + 4K \ ^ + 2£,tf J ^ + ' j 

x (2Q+q+q 1 q 2 - q+q 1 - g+g 2 ) + (off - diag.) , (bh) 
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where 



2\fi2 2A^w 1 6X2^2 



+ 71^ .9 +3 7^ + 



'2 



6 2 Q 2 



1 2 fii \ / 1 <pltt 2 



+2?hr+ 7777^" 77" + 



fii / V^2 4:K£ w u 



APPENDIX C 



In the present Appendix the expression of the operator iS d as well as the terms of the 
transformed expansion of the Hamiltonian (|2*4^) are given. 

In terms of RPA operators, the requirement (|25|) to eliminate the off-diagonal terms of 
order l/\/~N reads 

[H 2 ,iS] d = Y,tt„[QtQ„,iS] d = -(H 3 ) d . (CI) 



The solution of Eq. (jClj) can be found in the following form: 

iS d = iSY + isr (C2) 

with 

iSr = Y,\ v {Q+-QvMm-l) d , 

V 

iSr = ^* v (Q+-Q v ? d + p v (Q+-Q v ) 2 d (Q+-Q v ,)< l 

V 

+iAQt - QuUQt + Qu)l + 6„{Q+ - Qu'UQt + Q»)d 
+Pu(Qt - QvMQt + Q«)d(Qt> ~ Q»>)d - h -c- 

Inserting these expressions into Eq. ()C1|) we immediately find 

2 <j) v R v 2A vl — A u3 A v \ . , 

K = — 7 ==— ^- , a u = — , 7^ = -— . [LS) 
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The other three coefficients in Eq. (jC3|) are the solutions of a linear system of equations: 



Oj, 
Vl u r £l u 



\ 



/ 4 \ 

^%2 



A„ 5 



which gives 



Pu 



2A u2 n 2 u - A v4 (2nl - go - A^giQg 

O,,(4^-^0 
A v <£l v i — 2Q v (A v2 + A u4 ) 



With the condition (jCl|) one can eliminate the double commutators from Eq. 
obtain the following expression for the Hamiltonian H: 



(C4) 
to 



H = NH + (H 2 ) d + ^{[{HT + H: 
+2{H™ + HI* + H™)) +• 

J d 



(C5) 



It contains all corrections to RPA which are of the order 1/N. These corrections can be 
calculated by evaluating the expectation value of H. Of course, one could perform an 
additional unitary transformation to remove the remaining off-diagonal terms which are still 
present to order 1/N. This is, however, only required if one is interested in energies to 
order 1/N 2 . Since those are not needed for the present purpose all off-diagonal terms can 
be disregarded. 

Explicitly, the commutators appearing in Eq. (jC5|) are given by 



+ 



n 2 . 



[m—l 



a ■ 



[Hr,iS° 
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| Rv4>u{?>A vl - A u3 ) + R v ,$ v >{A v2 - A u4 ) 



n 



3/2 



x (2QtQ u + l) d (m - l) d + (off - diag.] 



3/2 
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[H™,iS™] d = 16j2{^ui(lu-a u ) + A u2 (35 u -P u ) + 3A u3 (3a u + lv ) 

+A vi {3(5 v - 8 V ) + A u5Pl/ yQ+Q+Q u Q u - 2Q+Q v ) d 
+ 16{(/? 2 - <J 2 )(A 13 - 3A n ) + (72 - 3a 2 )(A 12 - A 14 ) 
+2 Pl (A 12 + A u ) + 2A 15 (/3 1 + 50 + (ft - 50(^3 - 3A 21 ) 
+(71 - 3ai)(A 22 - A 24 ) + 2p 2 (A 22 + A 24 ) + 2A 25 ({3 2 + 5 2 )} 
x {2QtQtQ 2 Q l - QtQi - QtQ2) d + (off - diag.) . 

With these commutators we can write the next-to-leading order corrections in the following 
form: 

2 ^7T7T— »7T7T ( ^ ^ CXjCXi ^ ^ CKj + 1 ) , 

i([//3-,i^] d + [// 3 -,i^] d ) + # 4 7 

= ^K CT ,^(2g+g, + i) d (m-i) d 

V 

+ (off - diag.) , 

i[// 3 -,i^] d + // 4 7 = l -Y. v ^^^tQtQuQu-2QtQ,) d 

V 

(2g+g 2 + g 2 Q 1 -^Q+Q 1 ,) d 

V 

+ (off - diag.) , 

where 

8gg . 



VwiT— >7T7T 



^TTCT,,— >7T<7„ 
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6g£t 6^-^)(8^ + 2^-^j) 



x {nj(17^ - lOfi 2 ,) + {El - 20^)(20£ 2 + 7Q 2 ,)} 
2^-4^) 2 (4^-^,)(^-^) . 



X 



3Q?Q2(4Q2 _ fi2 / )(fi2 _ Q2)2 

3ft?Q|(ft 2 _ ^2)2 

{2(4^Q2, _ 5fi 2 fi 2 _ 1Q£ .2 Q 2 + 10 ^ 4) 



+ (4g "ll (4g j;^r (C6) 



4^2 - flj, 

Thus there are three sources of the residual interaction: pion-pion, pion-sigma and sigma- 
sigma couplings. Furthermore, we can evaluate the higher-order corrections to the pion mass 
as 

' u=l,2 

and to the sigma-boson mass as 
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